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■ The analysis of event shape variables enable studies of QCD and determinations of a s , 
but require predictions to next-to-leading order. The predictions of two next-to-leading 

Q_i| order programs are compared for current jet thrust, current jet-broadening, C parameter 

O ■ and jet-mass in deep inelastic scattering. 
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1. Introduction 



The large number of deep inelastic scattering (DIS) events collected at HERA permit the 
study of hadronic final states to high precision. Event shape variables describe the shape of 
jets of hadrons that are produced by the fragmenting quark. These have been investigated in 
e + e~ experiments, and can be used to extract the strong coupling constant independent of a 
particular jet-finding algorithm. Interest has recently been revived by the calculation of power 
corrections A similar analysis can be performed in DIS and compared to the e + e~ data. 

A number of general NLO programs have recently become available. The work presented 
here compares two programs, DISENTp] and DISASTER++0, interfaced through a common 
control library Q]. Both programs use the subtraction method^ and allow arbitrary (2+1) jet 
observables to be evaluated in the Breit frame of reference. The two programs were compared 
for both differential and mean distributions over a wide range in x and Q 2 to match with an 
existing ZEUS analysis]^]. 



2. Event Shape Variable Definition 

The event shape variables considered were current jet thrust T z f^\ current jet broadening B c , 
invariant jet mass pe and the C parameter. The variables are defined so as to be infrared and 
collinear safe to avoid divergences in the Monte Carlo integration and are calculated at parton 
level. 



Ei \n-Pi\ _ 


EiPz 


Ei \Pi\ 


Ei \p%\ 


Ei \n x pi\ 


EiP± 


2£* \Pi\ 


2£* \pi 


M 2 




(2E tot ) 2 ~ 





(1) 



PE ~ (2E tot ) 2 ~ 4(E, Eif 

C = 3(A!A 2 + A 2 A 3 + AiA 3 ) (4) 

where E tot is the total energy of all the particles in the current region. Ai 2,3 are the eigenvalues 
of the linearised momentum tensor, Q a >P, whose components are given by: 

QQ./3 _ EifijPi /\Pi\ /gs 

Ei\Pi\ 

The sum i is over all the partons in the current region of the Breit frame. The particles' three- 
momenta are given by pi, and their energies by E^. The momentum four- vector of the particles 
is denoted pf , and n is a unit vector aligned along the virtual boson axis (n = [0, 0, — 1]). 

In the Born approximation corresponding to the quark parton model (QPM), all the above 
variables are equal to zero, except thrust which is equal to unity. For consistency with other 
variables, it is convenient to define r z = 1 — T z , which goes to zero in the Born approximation. 

1 It should be noted that this is different from the standard definition of thrust T where the thrust axis has to 
be searched for by maximising the thrust in the expression. For T z , the thrust axis is defined by the virtual 
boson axis. 
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3. Program Comparison 

For both programs, DISENT and DISASTER++, events were generated in 17 kinematic bins of 
x and Q 2 as shown in Table [I]. The matrix elements produced by the programs are convoluted 

Table 1 

Kinematic bins used in the comparison. 
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1 


0.0006 - 


0.0012 


10.0 


- 20.0 


9 


0.0024 - 


0.0100 


160.0 


- 320.0 


2 


0.0012 - 


0.0024 


10.0 


- 20.0 


10 


0.0100 - 


0.0500 


160.0 


- 320.0 


3 


0.0012 - 


0.0024 


20.0 


- 40.0 


11 


0.0100 - 


0.0500 


320.0 


- 640.0 


4 


0.0024 - 


0.0100 


20.0 


- 40.0 


12 


0.0100 - 


0.0500 


640.0 - 


- 1280.0 
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0.0012 - 


0.0024 


40.0 


- 80.0 


13 


0.0250 - 


0.1500 


1280.0 


- 2560.0 
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0.0024 - 


0.0100 


40.0 


- 80.0 


14 


0.0500 - 


0.2500 


2560.0 


- 5120.0 
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0.0024 - 


0.0100 


80.0 • 


- 160.0 


15 


0.0600 - 


0.4000 


5120.0 - 


- 10240.0 
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0.0100 - 


0.0500 


80.0 


- 160.0 


16 


0.1000 - 


0.6000 


10240.0 


- 20480.0 












17 


0.2000 - 


0.8000 


20480.0 


- 40960.0 



with a parton density function (PDF) and the value of a s evaluated at the relevant scale. The 
CTEQ-4M [0 parton density function was used based on the MS renormalisation scheme. The 
renormalisation and factorisation scales, fiR and fip, were set to the momentum transfer Q. 

To ensure that the shape variables are infrared and collinear safe, an energy cut must be 
imposed so that the Monte Carlo integration does not diverge ||. It is required that at least 
3% of the total energy available to the current region, Q/2, has actually gone into the current 
region. 

Figures [I] to |] show the differential distributions for the four variables. The kinematic bins are 
arranged so that x increases going from left to right, and Q 2 increases going from bottom to 
top. The bin widths were chosen to match an existing ZEUS analysis. 

• For thrust in figure [l], a clear evolution in Q 2 is visible where as Q 2 increases, the high 
thrust population increases as the events become more collimated. Both programs agree 
very well over the whole kinematic range, and are consistent within statistical errors. In 
the high thrust region, where the events become very collimated, both programs encounter 
collinear singularities. In this region the calculations diverge (the highest thrust bin 
0.9 <T Z < 1.0 has a negative value indicating divergence). 

• The C parameter distributions in figure |2| show an increasing tendency to favour (2+1) 
jet-like events with increasing Q 2 . The expected fall off|5|] at C — | is reproduced by both 
programs. There is reasonable agreement between the two programs, particularly in the 
high Q 2 region. The same collinear singularities that cause thrust to diverge are visible 
as C — > 0; this appears either as a divergent point or a point with a large error. For 
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Figure 1. T z differential distribution. 
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Figure 2. C parameter differential distribution. 
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Figure 3. B c differential distribution. 
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Figure 4. p c differential distribution. 
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DISASTER++ a few kinematic bins show a point in the middle of the distribution with 
a large error indicative of some divergent behaviour; this is currently being investigated 
by the author of DISASTER++. 

• Jet-broadening in figure ^| shows the expected behaviour; the broadness of the jet reduces 
as Q 2 increases. The programs agree well in the high Q 2 region, however below a Q 2 of 
100 GeV 2 , there is some disagreement. 

• Figure f| shows good agreement for the jet-mass in most of the kinematic bins, although 
some of the medium Q 2 bins (160 — 320 GeV 2 ) by around 10% level. The errors in 
the low Q 2 regions are comparatively large indicating some problem in the numerical 
stability of the Monte Carlo integration in these regions for this observable. The numerical 
stability can, in general, be improved by increasing the number of points generated in 
each kinematic bin. 

Figure [| shows the mean distributions for the four variables as a function of mean Q, compared 
for the two programs. 

• Figure || a) shows the mean r z behaviour with increasing Q. This illustrates more clearly 
than the differential plots that the events become more collimated as < Q > increases. 
The different points at the same < Q > cover different x bins; a clear x dependence can 
be seen. Both programs show good agreement^] except for the x dependence. 

• The mean C parameter in figure [| b) shows the same behaviour as thrust with < Q >. 
The agreement is very good although there is a clear problem in the two lowest Q 2 bins. 
Both programs show minimal x dependence for this observable. 

• Figure |^ c) shows the mean jet-broadening. Above < Q > of 10 GeV the agreement is 
very good. Below < Q > of 10 GeV, DISENT shows a clear turnaround in the distribution 
which is not observed in DISASTER++. 

• The jet-mass mean distribution is shown in figure [|d). There is in general good agreement 
with both programs showing little x dependence. DISASTER++ has difficulties with 
divergences in two of the bins which require further investigation. 

4. Conclusions 

A comparison between two next-to-leading order programs, DISENT and DISASTER++, has 
been presented for four event shape variables in deep inelastic scattering. Differential and mean 
distributions have been considered. Both programs show good agreement over a wide kinematic 
range in the differential distributions, although jet-broadening shows some difference in the 
lower Q 2 region. The expected behaviour with increasing Q 2 is reproduced by both programs. 
The mean distributions show reasonable agreement except for jet-broadening where DISENT 
shows a clear turnaround in the distribution, which is not observed with DISASTER++. 

2 There was an error in the my code which rendered the errors in < r 2 > meaningless; they have been omitted 
for clarity. 
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Figure 5. Mean event shape distributions vs. < Q >. 
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